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ABSTRACT 

Aims. We aim to systematically investigate the cores of a sample of fossil galaxy groups and clusters (“fossil systems”) 
using Chandra data to see what hints they can offer about the properties of the intracluster medium in these particular 
objects. 

Methods. We chose a sample of 17 fossil systems from literature with archival Chandra data and determined the cool- 
core fraction for fossils via three observable diagnostics, namely the central cooling time, cuspiness, and concentration 
parameter. We quantified the dynamical state of the fossils by the X-ray peak/brightest cluster galaxy (BCG) separation, 
and the X-ray peak/emission weighted centre separation. We also investigated the X-ray emission coincident with 
the brightest cluster galaxy (BCG) to detect the presence of potential thermal coronae. A deprojection analysis was 
performed for fossils with s < 0.05 to resolve subtle temperature structures, and to obtain the cooling time and entropy 
profiles. We also investigated the Lx — T relation for fossils from the 400d catalogue to test if the scaling relation 
deviates from what is typically observed for other groups. 

Results. Most fossils are identihed as cool-core objects via at least two cool-core diagnostics with the population of weak 
cool-core fossils being the highest. All fossils have their dominant elliptical galaxy within 50 kpc of the X-ray peak, and 
most also have the emission weighted centre within that distance. We do not see clear indications of a X-ray corona 
associated with the BCG unlike that has been observed for some other clusters. Fossils lack universal temperature 
profiles, with some low-temperature objects generally not showing features that are expected for ostensibly relaxed 
objects with a cool-core. The entropy profiles of the z < 0.05 fossil systems can be described well by a power law 
with shallower indices than what is expected for pure gravitational processes. Finally, the fossils Lx — T relation shows 
indications of an elevated normalisation with respect to other groups, which seems to persist even after factoring in 
selection effects. 

Conclusions. We interpret these results within the context of the formation and evolution of fossils, and speculate that 
non-gravitational heating, AGN feedback in particular, could have had an impact on the ICM properties of these 
systems. 
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1. Introduction 

Clusters of galaxies are the largest gravitationally bound 
structures in the Universe, with an aggregate mass between 
a few times 10 13 M 0 and 10 15 M 0 . Consisting of galaxies, 
hot X-ray emitting gas (the intracluster medium-ICM) and 
dominated by dark matter, these objects are unique in their 
application to study both cosmology and astrophysics. On 
the lower mass end, these systems are sometimes called as 
groups, to indicate a smaller aggregation of galaxies. The 
distinction between groups and clusters is quite loose, with 
definitions in literature based on opti cal richness or simply 
a mass/temperature cut (e.g. IStott et alJl2012i) . 

Within the existing division of groups and clusters, lie 
an ostensibly ‘special’ class of systems called as fossil sys¬ 
tems; systems which are dominated by a single, bright, el¬ 
liptical galaxy and deficient in othe r br ight galaxies. The 
first fossil group was reported bv lPonman et all ( 1994 ) who 
showed the existence of a system at redshift 0.171 which 
had a large, X-ray halo associated with an elliptical galaxy. 
The choice of nomenclature ‘fossil’ is used to indicate a 
class of systems which is a remnant of galaxy merging, 
with hot X-ray gas and an elliptical galaxy the only rem¬ 
nants of this process. A ‘formal’ definition for a fossil sys¬ 


tem was provided bv I Jones et al.l (120031) who defined it as 
a spatially extended X-ray source with luminosity larger 
than 10 erg/s with an optical counterpart wherein the 
magnitude difference between the first and second bright¬ 
est galaxy is greater than or e qual to 2 in the Johnson 
R band. Simulations by iDariush et al.l (120071) show that 
systems selected by this criterion represent a class of ob¬ 
jects which assemble a higher fraction of their mass at high 
redshift as compared to non-fossils, further justifying their 
tag of being early-forming ‘old’ systems. This formal defi¬ 
nition, however , is open to contention as pointed out again 
bv IDariush et al.l (I2010D who argue via simulations that a 
magnitude difference of 2.5 between the first and fourth 
brightest galaxy identifies 50% more early formed systems. 

The most widely acknowledged formation scenario for 
these systems is one wherein dynamical friction causes 
galaxies close to the centre of the group to merge and 
leave behi nd a large elliptical g alaxy and a hot X-ray 
halo (e.g. iD’Onghia et al.l 120051) . An outstanding ques¬ 
tion is if the properties of fossil systems make them a 
special class of systems. Also, it is still unclear whether 
they are the final stages of mass assembly or are merely 
a temporary state in the formation of larger clusters. 
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lLa Barbera et al.1 (2009) argue that fossils are not special, 
and merely represent the final stages of mass assembly in 
a regio n witho ut enou gh surrounding matter. Simulations 
bv Ivon Benda-Beckmann et al. ' (120081 ) claim however, that 
the ‘fossil’ phase is a temporary one and the magnitude gap 
would end by renewed infall fro m the sur rou ndings . Adding 
to the confusion are results bv lProctor et al : (|201l h who ar¬ 
gue that the high mass to light ratios of these systems and 
their low richness are not predi c ted b y any current model 
of fossil formation. iPierini et all (l201lh argue that compact 
groups are the progenitors of low-mass fossils, while more 
massive fossil systems are formed by an early infall of mas¬ 
sive satellites, hinting that a singular explanation for the 
formation of fossils might be inadequate. 

Very few fossil systems have been thoroughly stud¬ 
ied, which adds to the lack of sound conclusions on their 
formation and evo l ution ( e.g. Mulchaev fc Zahhidofflll999t 
iKhosroshahi et all 12004 ISun et al.l l2004f) . X-rav sample 


tional details, are even rarer (e.g. 

Khosroshahi et al.l 20071 

Miller et al. 2012!; Harrison et al. 

2012j) and have focused 


more on global properties such as X-ray luminosity, tem¬ 
perature and mass. Till date, there has never been a de¬ 
tailed and focused investigation of the cores of a sample of 
fossil systems in X-rays, particularly with respect to their 
cool-core properties. 

The cores of galaxy clusters and groups offer impor¬ 
tant hints on the dynamical state of the ICM in the cen¬ 
tral regions, e.g. systems which are morphologically dis¬ 
turbed, generally have long central cooling times (CCTs) 
and elevated central entropies as compared to relaxed sys¬ 
tems which show much shorter CCTs and lower central en¬ 
tropies. The clusters with the strongest cooling, namely the 
strong cool core (SCC) clusters, show very short CCTs (< 
1 Gyr), low entro p ies (< 3 0 keV cm * 2 ), and central tempera¬ 
ture drops ()Hudson et al. 2010). This is not necessarily true 
on the group regime, where there are examples of systems 
which do not show any indication of cool g as in the cores de¬ 
spite short CCTs and low entropies (e.g. iBharadwai et al.1 
I2014H . Feedback also plays an important role in clus¬ 
ter/group cores as without it, there would be catas¬ 
trophic cooling with very high m ass deposition rates , 
in consist ent with observations 

Tamura et al.l J2001; 


Kaastra et al 


Peterson et al.,! ;2001fc 


;e.g. retei 
foJ 120011 


Peterson et al.1 


2003: ISanders et ah 120081) . Where do fossil systems fit in 


this picture? What hints can the cool core properties of 
these systems offer in terms of their formation and evolu¬ 
tion? These are the questions we attempt to address in this 
paper. 

This paper is organised as follows: Section [2] deals with 
the data and analysis. We present our results and discuss 
them in section [3] A short summary is presented in sec¬ 
tion [I] Throughout this work, we assume a ACDM cos¬ 
mology with = 0.30, 12 a = 0.70 and h = 0.70 where 
Hq = 100 h km/s/Mpc, unless stated otherwise. All errors 
are quoted at the 68% level. Log is always base 10 here. 


2. Data and Analysis 

2.1. Sample 

For this study, we selected a sample of 17 fossil systems and 
candidates spanning a large redslrift range (from ~ 0.01 to 
~ 0.4) with good quality archival Chandra data (f exp > 10 


ks) as required for exploring the properties under consider¬ 
ation here. Chandra was preferred over XMM-Newton due 
to its superior spatial resolution which makes it easier to 
probe the central regions which is the major focus of the 
work. Though the systems were categorized as fossil sys¬ 
tems through different techniques, qualitatively the basic 
selection criteria of a lower X-ray luminosity threshold and 
a magnitude difference between the first and second bright¬ 
est galaxies, is maintained for all of them. We also endeav¬ 
oured to study the Lx — T relation for those fossil systems 
selected from the 400d cluster sample (Sec. 13.61) for which 
we analysed XMM-Newton data for one more fossil group, 
which was not included for the other studies. Tabled] sum¬ 
marises the basic information of the sample with redshifts 
and literature sources. 


2.2. Basic data reduction 

Most of t he data reduction steps wer e similar to that 
follow ed in IBharadwai et al.1 (12015H and iBharadwajet al.l 
(1201411 and we describe them here briefly. Cl AO 4.43 with 
CALDB 4.5.0 was used for the data reduction for all but 
one object-cll416p2315. While the manuscript was under 
review, a very recent (2014), and much deeper observation 
of cll416p2315 was made available in the Chandra archives 
which we analysed with CIAO 4.7 and CALDB 4.6.7, to ac¬ 
count for changes in the ACIS time-dependent gain and to 
use the proper blank-sky background Hies. This is unlikely 
to have an effect on other objects in our sample which have 
much older observations, and is already accounted for in 
CALDB 4.5.0, and therefore we decided not to re-analyse 
them. Other than these changes, all the steps that were 
carried out were exactly the same. 

The chandra_repro task reprocessed the data and re¬ 
moved afterglows, created the bad-pixel table and applied 
the calibration. The lc_clean algorithm was used to fil¬ 
ter the soft proton flares, with the lightcurves visually in¬ 
spected to check for any residual flaring. Point sources 
were detected and excluded from further analysis with 
the wavdetect wavelet algorithm. Since our focus is on 
the bright, central regions, background subtraction for the 
spectral and surface brightness analysis was done through¬ 
out using the blank sky background files. The X-ray emis¬ 
sion peak (EP) was determined in an exposure corrected, 
point source subtracted image using the CIAO tools built 
into ds9 with a pixel scale of 2. 

For extracting spectra, we created annular regions with 
a minimum source counts threshold of 1000. This ensured 
that the annuli are sufficiently small and at the same time 
does not sacrifice too much on the signal to noise. The 
resultant spectra were fit with an absorbed APEC model 
in the 0.7 7.0 keV energy range, with the «h value taken 
from the built in «h ftools which outputs the value from 
the L eiden/Argentine/Bonn (LAB) survey dKalberla et al.l 
120051 1. The abundance was kept thawed for the spectral 
fit and in cases where we could not get a reasonable con¬ 
straint on the value, we froze it to 0.30. We used the 
lAnders fc Grevesse (Il98 9) abundance table throughout. 

The surface brightness profile (SBP) was obtained by 
centring on the EP in the 0.5-2.0 keV energy band. 
This SBP was then fit by a single or double beta model 


1 http://cxc.harvard.edu/ciao 
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Table 1. The sample of fossil systems. The starred entry represents the system which was only considered in the Lx — T 
study. Columns are (1) name of the fossil system, (2) redshift (3) literature sources. 


Name 

redshift 

Literature 

NGC 6482 

0.0131 

Khosroshahi et al. (2004) 

NGC 1132 

0.0232 

Mulchaev fe Zabludoff (1999) /Yoshioka et al. (2004) 

RX J0454.8-1806 

0.0314 

Yoshioka et al. 12004) 

ESO 306-G 017 

0.0358 

Sun et al. (2004) 

UGC 842 

0.045 

Voevodkin et al. 12008)/Lopes de Oliveira et al. (2010) 

RX J1331.5+1108 

0.081 

Jones et al. 12003) 

clll59p5531 

0.081 

Vikhlinin et al. 11999)/Voevodkin et al. (2010) 

cl2220m5228 

0.102 

Voevodkin et al. (2010) 

cll038p4146* 

0.125 

Voevodkin et al. (2010) 

cll416p2315 

0.138 

Jones et al. (2003)/Voevodkin et al. (2010) 

cl0245p0936 

0.147 

Voevodkin et al. (2010) 

cll340p4017 

0.171 

Ponman et al. (1994) /Voevodkin et al. (2010) 

RX J2247.4+0337 

0.199 

Vikhlinin et al. 11999)/Voevodkin et al. (2010) 

RX J0825.9+0415 

0.225 

Eigenthaler & Zeilinger (2009) 

RX J1256.0+2556 

0.232 

Jones et al. (2003) 

RX J0801+3603 

0.287 

Eigenthaler & Zeilinger 12009) 

RX J1115.9+0130 

0.352 

Eigenthaler & Zeilinger (2009) 

RX J0159.8-0850 

0.405 

Eigenthaler & Zeilinger (2009) 


(e.g. ICavaliere & Fusco-Femianoj['l97(f l given by: 

-3/3+1/2 

S = E 0 1 + ( — 


(1) 


£ = £, 


01 


1 + 


r-| -3/3i + l/2 


+ Eq2 


1 + 


-302 + 1/2 


( 2 ) 

where x Ci is the core radius. This in turn gave us the elec¬ 
tron density profile for the single and double beta cases 
as: 


n 


n 



( 3 ) 



1/2 

( 4 ) 


where r c is the physical core radius. The central electron 
density no for the single beta case is given by: 


n o = 


10 14 47rD A L> L CfVy 
El ) 


( 5 ) 


Here, N is the normalisation of the APEC model in the 
innermost annulus, £ is the ratio of electrons to protons 
(~ 1.2), Da is the angular diameter distance and L>l is the 
luminosity distance. El is defined as: 


El = 27 t / / x ( 1 + 

J — ooJ 0 


x 2 + L 


-30 


dxdl 


( 6 ) 


where R is the radius of the innermost annulus. 

_Similarly, for the double beta case, no (Appendix A of 

[Hudson et al.ll2010h is given by: 


no = 


10 14 47t(£ 12 LI 2 + Lh)D A D h CN 


£i 2 LI 2 EIi + LIiEI 2 


( 7 ) 


with the same definitions as before. LI; is the line emis¬ 
sion measure and is defined as: 


LI, 




( 8 ) 


To estimate the uncertainty on no, we conducted Monte 
Carlo simulations where the surface brightnesses were var¬ 
ied within their errors to generate new SBPs. These new 
SBPs were then fit again to a single or double beta model 
which was then used to determine no (the APEC normalisa¬ 
tion is also varied each time within the errors). The process 
was repeated 500 times to get a distribution for the values 
of no, and the standard deviation of the distribution gives 
the uncertainty on the measured no- 


2.3. Cool-core analysis 


The most robust paramete r to identify th e coo l-core nature 
of an object is the CCT dHudson et al.ll2010i) . a quantity 
that is dependent on the central temperature and the cen¬ 
tral density. But, considering that most (59%) of our ob¬ 
jects have a redshift >0.1, the determination of the CCT 
becomes problematic given that it is difficult to resolve the 
temperature profile for these objects in the central regions. 
Moreover, in some cases the data quality is not good enough 
to determine a temperature profile. Thus, in addition to the 
CCT, we used two other diagnostics, namely, Cuspiness a 
and concentration parameter csb, both of which have been 
sho wn to have a strong correlation with t he cooling time 
(e.g. lSantos et al.ll2010UHudson et al .1120101) . Indeed, recent 
studies of high redshift systems have used these t wo quan¬ 
tities to identify their C C nature (e.g. ISemler et al.i r2Q12: 
iPascut fe Ponman] 120151) . In our study, if a fossil satisfied 
two out of three CC diagnostics, it was classified as a CC 
system. 

The CCT is given as: 


CCT = t„„,(0) = 5 C +. a + nM 


2 ’ nt n A(To,Zo) 


( 9 ) 


where n;o and n e o are the central ion and electron densities, 
respectively, and Tq is the temperature in the innermost an- 
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nulus. As with iHudson et all (l2010h and iBharadwai et al.1 
(|2014fl we took the value of central density values to be the 
value at r = O.CMMrsoo- To estimate rsoo, we first determined 
the virial temperature (T v j r ) in a single region centred on 
the EP, which extended to the outermost boundary of the 
outermost temperature bin. The size of this region effec¬ 
tively corresponds to 0.3-0.5r 5O o for most objects, where 
the temperature gradient is not too significant (see Fig. [2] 
right, a ppendix IA1). We then used the following scaling rela¬ 
ti on by Evrard et al.l (1 19!) fill to determine r^nn as wa s done 
in lHudson et al.h2010f) aiid lBharadwai et alT (|2Q14h : 


^500 = 2 x 


Tvir V 

10 keV ) 


Mpc 


( 10 ) 


once again determined by Monte Carlo simulations as with 
n 0 - 

The concentration parameter csb (jSantos et al.ll2008lf is 
given as: 

E(r < 40 kpc) 


csb = 


E(r < 400 kpc) 


(17) 


where E is the integrated surface brightness within the de¬ 
fined radii (40 and 400 kpc). A CC system has csb > 0.075. 

As the z < 0.05 fossils are too close for estimating the sur¬ 
face brightness within 400 kpc, we did not determine csb 
for these systems. 


3. Results and Discussion 


To determine T v i r for C11038p4146 we used the XMM- 
Newton data sets which were analysed with the stan - 
dard pipelines as illustrated in lSchellenberger et al.l ( 2015 1. 
Spectrum was once again extracted in a single region 
centred on the EP. Note that temperatures determined 
from Chand ra and XMM -Newton can differ significantly 
as shown by ISchellenberger et al. l (1201511 and hence to re¬ 
main consistent, we used the formula given in that study 
to convert the XMM-Newton temperature to Chandra tem¬ 
perature: 

^°§10 (TxMm) = 0.889 X log 10 (Tchandra) (H) 

There are five fossil systems for which we could not ob¬ 
tain a temperature profile. Hence, for these objects we as¬ 
sumed T v ir to be the central temperature while constrain¬ 
ing the CCT. To estimate the uncertainties on the CCT for 
these objects, we conservatively assumed the lower bound 
of the temperature to be 0.6 times of T v j r (the largest cen¬ 
tral temperature decrement observed for nearby fossils, i.e. 
ES0306) and the upper bound to be 1.25 times T v ; r (the 
largest central temperature increment, i.e. NGC 6482). A 

CC system has a CC T < 7.7 Gyr. _ 

The Cuspiness a (|Vikhlinin et al.i r20071 is defined as: 


d log(n) 
d log(r) 


( 12 ) 


which in model parameters (IHudson et al.1 l2010h for the 
single beta model is expressed as: 


a = 



(13) 


and for the double beta model is: 


q = 3r 2 S i2 LI 2^H Cl 2 ^i + LIiff 2 r C2 2 &2 
S12LI261 -f- LIi6 2 


(14) 


Here, Ei 2 ,LIi and LI 2 are as defined for the CCT, bi and 
b i are defined as: 



and 



(16) 


respectively. In both single and double beta cases, r = 
0.04f5oo- A CC system has a > 0.5. The error on a is 


3.1. Cool core properties 

Table [2] summarises the CC properties of the fossil systems. 
We find that 14 out of 17 systems, i.e., 82% are cool core 
systems, as evidenced by at least two out of the three CC 
diagnostics. If one uses the CCTs, we observe a range of 
values from a minimum of 0.124 Gyr to a maximum of 
24.2 Gyr. The sub-classifications of the fossil systems using 
just the CCTs are 29% SCC (CCT < 1 Gyr), 53% WCC 
(1 < CCT < 7.7 Gyr), and 18% NCC (CCT > 7.7 Gyr). 

As we have already pointed out, it is not always pos¬ 
sible to determine the CCT and this becomes especially 
true for X-ray survey data (e.g. eROSITA). Thus, we also 
state the effect of using one CC parameter alone. With only 
the cuspiness, we identify 76% of the fossil systems as CC 
systems. If one uses the concentration parameter alone, we 
observe that 67% of the systems are selected as CC systems 
(though it’s quite likely that the 2 < 0.05 fossils will also 
have csb > 0.075, in which case the fraction would go up 
to 76%). 

Thus, if one were to simply use a single quantity to 
identify CC groups, conservatively, one would get a lower 
limit on the fossil CC fraction of 67% and an upper limit of 
82%, i.e. qualitatively there are indications that a majority 
of fossil systems are CC systems, albeit as the cooling times 
indicate, most seem to be WCC systems. Fig. [T] shows the 
relation between the CCT-a and CCT-csb and despite the 
scatter, it is clear that the anti-correlation trend between 
the quantities is quite strong with Spearman correlation 
coefficients of -0.91 and -0.80 respectively. 

Assuming no significant effect of selection bias on the 
sample, to the best of our knowledge, this is the first 
demonstrated result that suggests that most fossils are 
in fact cool-core systems based on observable diagnos¬ 
tics. The fraction of cool-core fossil systems is marginally 


d (77%. Bharac 

wajet al. 20141) and cluster sam- 

Hudson et al. 

2010), though we do point out 


that studies have shown that cool-core fractions are gen¬ 
erally over : estimated_Jn flux-limited samples of clusters 
(e.g. IHudson et ahl l201(ll : lEckert et al.l l201lt iMittal et al.l 
120111) . The high cool-core fossils fraction is not totally un¬ 
expected as it is generally assumed that fossils are among 
the most relaxed clusters/groups in the Universe and hence 
should have pronounced cool-cores with very short CCTs. 
What is more intriguing is the fact that the WCC fraction 
is the highest among the three classes, and the population 
of WCC and NCC fossils combined together (collectively 
called non-strong cool-cores; NSCCs) exceed the SCC fos¬ 
sils population. This is not what one would naively expect 
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Table 2. Cool-core diagnostics of the fossils. The columns are (1) name of system, (2) coordinates (J2000) of the EP, 
(3) virial temperature of fossil, (4) central cooling time, (5) cuspiness, (6) concentration parameter, (7) cool-core or not, 
via any two diagnostics. Starred entries represent fossils for which the CCT was determined using T v j r . 


Name 

EP (RA/DEC) 

Tvir (keV) 

CCT (in Gyr) 

a 

CSB 

CC or not? 

NGC 6482 

17 : 51 : 48.85 

+ 23 : 04 : 19.34 

rj 7^+0.0078 

u -' '-'-0.0079 

0 124 +0 ' 015 
u.iz ^_ 0 013 

1 46 +0 - 01 

^^-O.Ol 

- 

YES 

NGC 1132 

02 : 52 : 51.79 

- 01 : 16 : 29.13 

1 90 +O.OI 5 

l.ZO—o.016 

1 iq+0.18 

J-O—0.14 

0 Q 6 Q +0 066 

u.yuy_ 0 .066 

- 

YES 

RX J0454.8-1806 

04 : 54 : 52.34 

- 18 : 06 : 55.57 

2 Sfi + 0 07 

z.OU—0.07 

901 +1.38 
^• O1 -0.78 

PI QQQ+0-020 

U.OOO_o 020 

- 

YES 

ESO 306-G 017 

05 : 40 : 06.75 

- 40 : 50 : 10.80 

0 71 +0.07 
' 1 -0.05 

0 4Q7+0 - 064 
/ —0.057 

0.7901+22 

- 

YES 

UGC 842 

01 : 18 : 53.67 

- 01 : 00 : 06.75 

1 Q1+0.07 
± ' y - L -0.07 

1 7 ^+ 0-22 

1 - 'O-0.18 

1 . 02 ++ 

- 

YES 

RX J1331.5+1108 

13 : 31 : 29.65 

+ 11 : 07 : 57.44 

n 74 + 0-02 

u - '^-0.03 

1 Q 1 + 0- 38 
iU± -0.23 

1 qq+ 0.11 
l.OO-o.n 

0 11 n + 0 006 
U. 11 U _0 006 

YES 

clll59p5531 

11 : 59 : 52.20 

+ 55 : 32 : 06.29 

1 70 + 0 .O 6 

1 u - 0.02 

0-261™ 

1 70 +O. 2 O 
‘ *- 0.20 

0 285 -1 " 0,005 

U.ZOO_o 005 

YES 

cl2220m5228* 

22 : 20 : 08.64 

- 52 : 27 : 50.33 

4.12+;+ 

20.9+11 

0.264+;°+ 

0 D62 +0 006 
U.UOZ_o 006 

NO 

cll416p2315 

14 : 16 : 27.39 

+ 23 : 15 : 22.59 

3.98±°+ 

4-731+1 

0.5481+53 

0 05Q" 1 ” 0 002 
U.UOy_Q 002 

YES 

cl0245p0936* 

02 : 45 : 48.83 

+ 09 : 36 : 37.30 

2 4> : i + 0 - 38 

4.99+'+ 

0 7Qf ) +0 - 028 
u. < yO_0.028 

0 11 fi + 0 - 012 
U. 110 _Q 012 

YES 

cll340p4017* 

13 : 40 : 32.70 

+ 40 : 17 : 39.26 

1 42+0- 12 
i-^_o.io 

q 9 Q+ 0.53 

o.z,y _ 0 44 

0 620 +0 ' 025 
u.ozu _ 0 025 

0 17Q + 0 - 011 

0.1 /y_ 0 .oii 

YES 

RX J2247.4 + 0337* 

22 : 47 : 27.88 

+ 03 : 37 : 00.22 

2-43++ 

15 0+3-58 
AO.y_2.84 

0-3591++ 

0.060+;“* 

NO 

RX J0825.9+0415 

08 : 25 : 57.83 

+ 04 : 14 : 47.12 

4.481+9 

5 lfi + 0 - 76 

O.IO— 0.55 

0.460+;+ 

0 . 0861 ++ 

YES 

RX J1256.0 + 2556* 

12 : 56 : 02.25 

+ 25 : 56 : 35.38 

q 1 q+0.75 

O. 1 O — 0.51 

9 4 0+5.15 
Z^l.Z-s 84 

0 25f) +0 ' 045 

U.ZOU_o.o45 

0 n74+0 - 010 

u.u/4_ 0 010 

NO 

RX J0801+3603 

08 : 00 : 56.81 

+ 36 : 03 : 23.67 

7 40 +O .29 
'-^O _ 0 30 

9 9 C+O .32 
^•^O_ 0 .24 

n 099 + 0.039 
U.O^_o 039 

0 1 2 n +0 003 

u.izu_ 0 .oo3 

YES 

RX J1115.9+0130 

11 : 15 : 51.80 

+ 01 : 29 : 55.48 

7 -1 q+0.26 

1 • ^-0.26 

0.629++1 

1 20 +003 
T.z,u_0.03 

0 224" 1 " 0 ' 005 
U.ZZ41_o 005 

YES 

RX J0159.8-0850 

01 : 59 : 49.36 

- 08 : 49 : 59.45 

8-46+33 

0.598++° 

1 io + 002 

1-J-U-0.02 

H 9Qq+0.006 

u.z,z,o_ 0 .oo6 

YES 


from ostensibly relaxed objects such as these and are first 
indications of an impact of non-gravitational processes on 
the fossil ICM. 

3.2. EP-BCG/EP-EWC separation 

We also quantified the dynamical state of the fossil sys¬ 
tem via two properties, namely the emission peak-emission 
weighted centre (EWC) separation and the central elliptical 
galaxy-EP separation, i.e. the BCG-EP separation. These 
are reasonable indicators of whether the system is relaxed 
or not, as large values f or these n u mbers would indicate a 
disturbed system (e.s. iHudson et all 1201 Oft . 

The EWC was determined iteratively in a point source 
subtracted, exposure-c orrected i mag e by st arting off from 
the EP as was done in iHudson et ali (l2Q06h . The positions 
of the BCG were taken from literature and was confirmed 
from the NASA extragalactic database (NED0- 

Table [3] shows the EP -EWC and EP-BCG separation. 
All fossils have their BCG located within 50 kpc of the X- 
ray peak. The EP-EWC peak shows a wider range of values 
with four fossils (24%) having a separation > 50 kpc. There 
is an anti-correlation trend of the EP-EWC separation with 
the Cuspiness with a Spearman correlation coefficient of 
-0.86, i.e. cool-core fossils have a smaller separation and 
vice-versa (Fig. [2j left) 

Despite most fossils having an EP-EWC separation less 
than 50 kpc, again such a large variation in EP-EWC values 
is not naively expected for ostensibly relaxed systems such 
as fossils. 

3.3. Temperature Profiles 

Figure [2j right and appendix [Al shows the scaled temper¬ 
ature profiles for the fossil systems for which we could ob¬ 
tain one (12 objects). There’s a great degree of variation in 

J http://ned.ipac.caltech.edu/ 


the profile shapes with temperature drops and flat shapes 
for the fossil systems and we fail to pinpoint a “typical” 
temperature profile for the fossils. The high redsliift, high 
temperature systems show temperature profiles that are ex¬ 
pected for a SCC/WCC cluster, i.e. indications of a tem¬ 
perature drop for the SCC ob jects a nd a relati vely flat pro¬ 
file for the WCC objects (jHudson et al.ll2010h . NGC 6482 
shows a profile that rises all the way to the centre which has 
been reported thoroughly in previous studies of the object 
(iKhosroshahi et al]l2004f) . NGC 1132 shows an almost a flat 
profile with barely a decline towards the centre. ESO306, 
RX J0454.8-1806, and UGC00842 lacks a smooth tempera¬ 
ture decline towards the centre which while potentially in¬ 
dicative of a cool-core, shows that it’s a truncated one. One 
possibility for the shape of such a temperature profile could 
be emission from the BCG, which we discuss in Sec. 13.41 
A deprojection analysis ISec. 13.51) for the z < 0.05 systems 
indicates the imprint of a heating process such as AGN. 
Indeed, only clll59p5531 shows the “classic”, peaky tem¬ 
perature pro file observed for low-temperature grou ps with 
a SCC (e.g. ISun et ahll2009t lEckmiller et al.lDOl ill versus 
that for SCC clusters which show f lattening at intermedi¬ 
ate radii, fe.g. IVikhlinin et al.ll200fil) . 

3.4. Potential emission from the BCG 

As shown in Table [3] most fossils have a BCG which is 
very close to or on the X-ray peak. The source detection 
tool wavdetect manages to detect sources which seem to 
be coincident with the BCG in most cases and we explored 
whether part of the emission in the core regions comes from 
the BCG itself, i.e. if they are coronae which consists of gas 
of stellar origin. This could in principle also explain the 
truncated “cool-cores” seen in objects such as UGC00842, 
ESO306, and RX J0454.8-1806. To study this emission, 
we focused our attention only on the fossil systems with 
z < 0.1, as it becomes harder to resolve potential coro¬ 
nae for high redshift objects. Using the regions generated 
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Fig. 1. Plot of CCT vs. cuspiness (left) and CCT vs. concentration parameter (right). The horizontal lines represent the 
cut for being classified as a cool-core system. The vertical line represents CCT = 7.7 Gyr. 


Table 3. EP-BCG, EP-EWC separation. The columns are (1) name of system, (2) separation between BCG and EP in 
kpc, (3) separation between EP and EWC in kpc. The uncertainty on these quantities correspond to 1 which in turn 
corresponds to 0.27 kpc for the nearest fossil, and 5.41 kpc for the furthest one. 


Name 

BCG-EP (kpc) 

EP-EWC (kpc) 

NGC 6482 

0.212 

1.16 

NGC 1132 

0.785 

5.90 

RX J0454.8-1806 

3.22 

13.3 

ESO 306-G 017 

0.842 

13.3 

UGC 842 

0.960 

2.87 

RX J1331.5+1108 

1.71 

14.8 

clll59p5531 

3.77 

4.40 

cl2220m5228 

3.71 

62.8 

cll416p2315 

0.964 

31.6 

c!0245p0936 

2.78 

32.4 

c!1340p4017 

3.46 

30.4 

RX J2247.4+0337 

2.93 

59.2 

RX J0825.9+0415 

3.58 

52.0 

RX J1256.0+2556 

6.6 

61.2 

RX J0801+3603 

1.29 

29.5 

RX J1115.9+0130 

7.88 

21.4 

RX J0159.8-0850 

4.31 

19.5 




Fig. 2. Left:P\ot of EP-EWC separation and cuspiness. Right: Scaled temperature profiles colour-coded as a function of 
virial temperature. 
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by wavdetect we extracted spectra which were fit to an 
absorbed APEC model. The regions used for the spectral 
extraction are of the order of a few kpc, which is the typical 
size of a corona. We kept the temperature, abundance, and 
normalisation as the free parameters. The best-fit values 
were used to determine the X-ray luminosity using Xspec 
in the 0.5-2.0 keV band, which was then compared with the 
ROSAT luminosity of the fossil, to determine the strength 
of this “corona” emission. 


Our estimated temperatures from the spectral fit is in 
most cases much higher (> 1 keV) than that is typically 
observed for the corona class. Moreover, the X-ray luminos¬ 
ity of the emission is in the order of 10 41 erg/s which is an 
order of ma gnitude h igher than what is typically observed 
for coronae ( Sun et aTJl2007t ) and can account for as much 
as 17% of the total fossil ROSAT luminosity. Lastly, the 
abundances are much lower (median of 0.28) than w hat has 
been reported for true coronae 10.8. ISun et al.ll2007 ') which 
gives indications that the detected emission might not be 
completely of stellar origin, though we do point out that 
there is evide nce from li terature for low-metallicity coronae 
as well (e.g. 1 Sanders et al.llWl4bl) . One strong possibility 
is that the stellar content has mixed with the surrounding 
dense ICM. Table [4] shows the properties of the emission 
associated with the BCG. 


3.5. Deprojection analysis of z < 0.05 fossils 

In order to explore thermodynamic properties of fossil sys¬ 
tems in detail, we focused in particular on the low red- 
shift (z < 0.05) fossils (five in total, namely, NGC 6482, 
NGC 1132, ESO306, RXJ 0454.8-1806, UGC 00842) and 
performed a deproj ection analysis . To do this, we adopted 
the MBPROJ code (ISanders et" al.ll2014ah whi ch is based on 
t he su rface brightness deprojection method of lFabian et al.1 
( 1981 ) and is used to estimate the ICM thermodynamical 
properties from the surface brightness profiles. This method 
potentially has the advantage of resolving the extreme in¬ 
ner regions of the galaxy group which cannot be probed 
with traditional spectral analysis in the absence of very 
high quality data. Additionally, this is a good opportunity 
to test the code on low temperature groups to see how well 
the results agree with spectral analysis and test the usabil¬ 
ity of the code for situations where spectral analysis would 
not be possible (e.g. survey data). 

For purposes of brevity we refrain from describing the 
code in detail and guide the reader to the aforementioned 
paper. The relevant information for this particular study 
are as follows: 


— Surface brightness profiles were extracted in the 0.5-1, 
1-1.5, 1.5-2.5, 2.5-7 keV bands. 

— We focused only on the “null potential” case for this 
particular study. This particular case does not assume 
hydrostatic equilibrium and essentially corresponds to 
low resolution spectral fitting. The advantage of this 
case is that it prevents the thermodynamic profiles from 
being biased due to assumption of hydrostatic equilib¬ 
rium. The null potential case can be used to estimate 
thermodynamic properties with and without linear in¬ 
terpolation of the temperatures between the bins. The 
non-interpolation case assumes the temperature is the 
same in groups of three density bins. For the scientific 


discussion here, we focused only on the profiles which 

do not involve linear interpolation of the temperatures. 

— The abundance is assumed to be 0.3 Z Q throughout. 

Figure [3] shows a comparison of the temperature profiles 
from the spectral analysis and MBPROJ. The code does a 
good job of recovering the shape of the temperature profiles 
though some deviations are seen for NGC 1132 and NGC 
6482 which could be due to the extremely low temperatures 
of these systems (1.2 and 0.75 keV respectively), and possi¬ 
ble projection effects. For the hotter groups, the agreement 
is much better. Overall though, MBPROJ seems to be a 
competitive code even on the group regime and would be 
particularly useful for eROSITA data where most systems 
would be galaxy groups and would lack suffi cient co un ts fo r 
a direct spectral temperature estimate dBorm et al.ll2014| ). 

Figure [I] shows the cooling time and the entropy pro¬ 
files of the fossil systems. The cooling time profiles can be 
well-described by a powerlaw with an index varying from 
1.06±0.03 to 1.20±0.05 for the five systems. iHudson et al.l 
(120101 ) define the cooling radius as the radius at which 
f cool = 7.7 Gyr. For these five fossils, namely, NGC 1132, 
RX J0454.8-i806, ESO306, NGC 6482, and UGC 00842, 
this corresponds to 43 kpc, 36 kpc, 51 kpc, 54 kpc, and 
30 kpc respectively. However, as shown in Fig. [3l the re¬ 
gion with the cool gas (if present) is generally much smaller 
than the cooling radius. MBPROJ detects cool gas out to 
3 kpc for NGC 1132, 12 kpc for RX J0454.8-1806, 9 kpc for 
ESO306, 0 kpc for NGC 6482, and 19 kpc for UGC 00842. 
Thus, all these systems seem to lack a group-sized cool-core 
and the region of cool gas is not commensurate to the region 
of short cooling times. As argued above, it seems unlikely 
that the cool gas component is emission purely associated 
with the BCG. 

The derived entropy profiles can also be well-described 
by a powerlaw of the form K (r) = K mo x (r/100 kpc) m as 
has been used in studies such as e.g. lCavagnolo et all ((2009) 
and [Panagoulia et al. (120141 ). When the profile is fit with 
the above function, a powerlaw index that is substantially 
lower than 1.1 is obtained; which is what is expected fo r 
p ure gravitat ional processes (e.g. iTozzi fe Normanl I2001D . 
ICavagnolo et al.l (. 2009 ) also fit their entropy profiles with 
a powerlaw model with an additional component for a non¬ 
zero core entropy, i.e. K{y) = K 0 + K wo x (r/100 kpc) m . 
We also added this component to see if it has an effect on 
the fitted entropy profile. As shown in Table [5J the index 
is still much lower than 1.1, and in some cases, A'o takes 
unphysical n egative values. Qualitati vely, this agrees with 
the results of lPanagoulia et akl ( 2014 ) who showed that the 
entropy profiles of a sample of galaxy groups and clusters 
can be sufficiently described by only a powerlaw without 
indications of any core-flattening. We also compared the 
scaled entropy pr ofiles to the baseline entropy profile of 
IPratt etahl ( 2010 ). which is defined as: 

K(R)/K 500 = 1.42(i?/ A 500) 1 ' 1 (18) 

where K 500 is as defined in IPratt et all (|2010l ) and is de¬ 
pendent on M 5 oo- T o estimate Msm we used the M 50 0 — T 
scaling relation from lArnaud et alJ (2005). As is evident in 
Fig. S] there is a significant entropy excess over the base¬ 
line entropy profile for the fossils within the radial ranges 
probed here, and the profile is considerably flatter as com¬ 
pared to the baseline profile. The shallow slopes, and the 
elevated entropies over a substantially large radial range 
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Table 4. X-ray emission coincident with the BCG for z < 0.1 fossils. Columns are (1) Name of group, (2) Temperature, 
(3) metallicity, (4) X-ray luminosity from Xspec, (5) “Corona” luminosity as percentage of ROSAT luminosity 


Name 

T 

z e 

Lx(0.5 - 2.0keV)(10 41 erg/s) 

% of ROSAT group emission 

NGC 6482 

0.951°'Q 22 

0.28±°;°! 

1.05 

17% 

NGC 1132 

0.82lH 6 8 

O-llto.04 

1.78 

3.8% 

RX J0454.8-1806 

1 04 +0 ' 16 
1 - u ^-0.08 

0.40t°;g 

1.06 

1 .2% 

ESO 306-G 017 

i oq+ 0.04 
1.0U_o.o5 

0 2Q+ 0,12 

U.zy_ 0.09 

6.46 

2.4% 

UGC 842 

1 QQ+ 013 

x.oy_ 0 .07 

0 21 +0,n 
u - zl -0.07 

1.78 

7.1% 

RX J1331.5+1108 

Pi qq+0-063 

u.oy_o.o34 

0.15tg;* 4 

1.84 

2.9% 

clll59p5531 

i qq+0.02 

i.zO-o.02 

0.38t°;° 6 5 

17.2 

15% 


is indicative of a strong influence of non-gravitational pro¬ 
cesses on the properties of the ICM, which we discuss in 
greater detail in Sec. 13.71 


3.6. Ax — T relation for 400d fossil systems 

The Lx — T relation is an interesting scaling relation 
which can give insights into the baryonic physics at play 
in the ICM as the physics affects both quantities, T and 
Lx- For example, clusters with cool-cores/non-cool cores 
have Lx — T relations with different slopes and normali¬ 
sations, and excising the core regions result in a lowered 
intrinsic scatter (e.g. iPratt et ahl 120091 iMittal et al.l 1201 It 
iMaughan et al.ll2012ll . Note that along with cool-cores and 
possibly AGN feedback, selection effects can also play a sig- 


been shown bv several authors before (e.g. Ikebe et al.i2002: 

Stanek et al.l 

200(1 Pacaud et al. 

2007| 

Mantz et al. 201C; 

Mittal et al. 

20111 Lovisari et al. 

2015 

Bharadwai et al. 

201jl also see Giodini et al. 2013 

section 7, for a review). 


Thus, simply using a sample of clusters/groups to construct 
scaling relations without accounting for potential selection 
effects would not result in an accurate description of the 
underlying properties of the objects. In this work, there are 
seven fossil systems from the 400d catalogue (out of 12). 
The 400d catalogue has fairly well-defined selection criteria, 
making it relatively easy to account for selection effects, and 
study if the Lx — T scaling relation for fossils show features 
that are different to non-fossils. To compare the results 
here, we used the groups in the [Blraradwai et al.1 (|2015h 
study. To maintain consistency with that study, we used the 
best-fit virial temperatures and abundances of the fossils, 
to convert the catalogue ROSAT (0.5-2.0 keV band) lumi¬ 
nosities into bolometric luminosities (0.01-40 keV band). 

To determ ine the scaling relat ion, we used the BCES 
(Y|X) code bv lAkritas fc Bershadv flOOfii ! and the fits were 
performed in Log space using the fitting function: 


/ L x (0.01 -40keV) \ _ / T \ 

\0.5 x 10 44 erg s -1 J ° \3 keVy 


(19) 


This gives a slope of 2.09 ±0.24 and normalisation (log) 
of 0.29 ± 0.04. The observe d scaling relation for groups in 
the lBharadwai et~ahl (|2015li has a slope of 2.17 ± 0.26 and 
normalisation of —0.01 ± 0.09. To ensure that our choice of 
region for the temperature determination was not biasing 
our results, we performed a test by estimating the temper¬ 
atures in a fixed region of 3 arcmin for all the 400d fossils 
(effectively corresponding to a median of ~ 0.5rsoo) and re¬ 


did the scaling relation. No significant changes were found. 
The observed scaling relation gives us first indications that 
fossils seem to be more X-ray luminous for a given T (Fig.O 
left), relative to non-fossils, though this could simply be due 
to selection effects. 


To investigate this further, we performed an additional 
test where we froze the slopes for both the fossils and groups 
to 3.0, and only fit the normalisation for both samples. 
The value of 3.0 for the slope was chosen as generally most 
groups and clusters can be described well by this particular 
slope, after factoring in selection effects. The fossils now 
have a normalisation of 0.38 ± 0.11 and the groups sample 
have a normalisation of 0.22 ± 0.13. We now proceeded to 
remove the influence of selection effects on the normalisa¬ 
tion of the scaling relation for both samples by generat ing 
mock samples of objects like in iBharadwai et al. (l2015f) by 
varying the input normalisations. For each mock sample, 
the slope was always fixed to 3, and the intrinsic scatter 
fixed to the observed values. Flux and luminosity cuts were 
applied to both the 400d sample and the groups sample to 
match the true sample of objects. For the 400d fossil sam¬ 
ple, th e selection criteria were taken from IVoevodkin et al.l 
(2010), and can be described as follows: a lower flux cut 
of 1.4 x 10 -13 erg/s/cm 2 , an upper redshift cut of 0.2, 
and a lower luminosity cut of 10 43 erg/s. Here, fluxes and 
luminosities are in the ROSAT (0.5-2.0 keV) band. After 
performing the bias-corrections, the fossils scaling relation 
has a normalisation of 0.30 ± 0.10 vs. 0.0078 ± 0.13 for the 
groups, indicating that the large normalisations seem to 
persist even after accounting for selection effects. Despite 
a nearly 2.3er higher normalisation for fossils with respect 
to non-fossils, we still treat this finding only as an indi¬ 
cation, because on top of the statistical uncertainties on 
both quantities (which is quite large currently for the 400d 
fossils), there could also be an additional systematic un¬ 
certainty introduced due to differences in the luminosity 
determination in the different parent catalogues. Secondly, 
the archival nature of this sub-sample of 400d objects is 
biased towards systems lacking a SCC (table [2]) and as¬ 
suming that the remaining systems are SCC, then adding 
these objects could potentially increase the normalisation 
for the scaling relation, and the difference between the fos¬ 
sils and the groups sample would be higher than what we 
demonstrate here. We plan to explore this in greater detail 
in a future study of fossils scaling relations. 
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Fig. 3. Temperature profiles of the fossil systems, from spectral analysis and MBPROJ. Clockwise from top left: NGC 
1132, RXJ 0454.8-1806, NGC 6482, ESO306, UGC 00842. The spectral temperature profiles were also determined for an 
abundance of 0.3 -Zq. 


3.7. Discussion 

Though our results concern mostly the cores of fossil sys¬ 
tems, these properties offer us interesting insights into their 
formation and evolution. Fossil systems not having experi¬ 
enced any recent major merger activity should have a pro¬ 
nounced cool-core with very short CCTs and ‘typical’ tem¬ 


perature profiles with a clear central temperature drop. The 
distribution of the CCTs and the shapes of the observed 
temperature profiles does not however seem to support this 
idea. Given that the BCG-EP separation is not large and 
that most fossils still have an EP-EWC separation less than 
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R/R500 


Fig. 4. Top: Cooling time profiles of z < 0.05 fossils. Middle: Entropy profiles of z < 0.05 fossils. Both quantities are 
obtained from the non-interpolated temperatures estim ated from the MB PROJ code. Bottom: Scaled entropy profiles. 
Dashed line represents the baseline entropy profile from lPratt et al.l (1 2010 *). 


50 kpc, it seems unlikely that invoking mergers can easily 
explain the thermodynamic properties of the fossils. 

A lot of these properties can be explained due to a 
non-gravitational process, in particular the AGN activity 
at the centres, which we speculate was much more power¬ 
ful at an earlier epoch than it is now. Simulation results by 


e.g. iDariush et al.l (|2007l) and [Dfaz-Gimenez et ahl (1201 lH 
point out that fossils assem ble most o f t heir masses at high 
redshifts. Interestingly, iDfaz-Gimenez et al.1 1 2008 1 point 
out that while the virial mass of the groups were assembled 
at high redshifts, much of the BCG mass is assem ble d at 
much later epochs. Recent results bv 1 Mendez-Abreu et al.l 
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Table 5. Best-fit results for Entropy profiles of z < 0.05 fossils, without and with core component. Columns are (1) 
Name of group, (2) Slope, (3) Normalisation, (4) core entropy 


Name 

m 

normalisation (Rioo,keV end) 

Ko keV end 

NGC 6482 

0.74 ±0.04 

103.3 ± 13.12 

— 

NGC 1132 

0.75 ±0.02 

138.9 ± 6.08 

— 

RX J0454.8-1806 

0.72 ±0.02 

184.2 ±9.83 

— 

ESO 306-G 017 

0.83 ±0.02 

170.1 ±9.00 

— 

UGC 842 

0.75 ±0.02 

224.49 ± 8.45 

— 

All 5 groups 

0.80 ±0.01 

146.1 ± 5.41 

— 

NGC 6482 

0.91 ±0.05 

139.7 ± 16.3 

1.24 ±0.29 

NGC 1132 

0.70 ±0.04 

138.9 ±5.96 

-1.38 ± 1.30 

RX J0454.8-1806 

0.72 ±0.05 

184.5 ± 11.1 

-0.14 ±2.38 

ESO 306-G 017 

0.90 ±0.04 

168.9 ±8.25 

1.83 ±0.84 

UGC 842 

0.70 ±0.03 

225.1 ± 7.72 

-2.43 ± 1.53 

All 5 groups 

0.84 ±0.02 

147.5 ± 5.38 

0.74 ± 0.36 


Table 6. Comparison of the “bias-corrected” scaling relations for the 400d fossils and the lBharadwai et al.l ( 2015 ) sample 
by freezing the slope to 3. 


Sample 

Slope 

Normalisation (observed) 

Normalisation (bias-corrected) 

400d fossils 

3 

0.38 ±0.10 

0.30 ±0.10 

Groups sample 

3 

0.22 ±0.13 

0.0078 ±0.13 




Fig. 5. Left: Observed Lx — T relation for 400d fossils and groups sample. Right: Comparison of 400d fossil Lx — T 
relation and groups relation with frozen slope and bias-corrected normalisations. 


(2012) conclude that the BCGs of fossils underwent ma¬ 
jor dissipational mergers at earlier epochs, albeit most of 
the stellar content was assembled via dissipationless ma¬ 
jor mergers at later times. Thus, we speculate that dur¬ 
ing the relatively ‘gas-rich’ phase of mass assembly, the 
SMBH was strongly fuelled which led to strong AGN feed¬ 
back such that an almost isothermal central region, or a 
centrally rising temperature profile was formed. When this 
phase ended, dissipationless major mergers continued to 
build the BCG, but there wasn’t enough gas to fuel the 
SMBH strongly, pushing it into a state of relative quies¬ 
cence. Eventually, for some objects, a small cool-core region 
was formed probably from the mixture of stellar mass ejecta 
and the surrounding ICM, not quite unlike thermal coro- 
nae seen in groups/clusters. This cool gas started feeding 
the SMBH, and is now probably involved in a self-regulated 
feedback cycle which leads to a much lower AGN activity. 
Four out of the five z < 0.05 fossils show extremely weak ra¬ 


dio sources (iBharadwai et alJl2014t iMiraghaei et al.ll2014t 
ICondon et al.lll998l) indicating that the AGN activity is def¬ 
initely not strong currently, giving credence to this_tlieory. 
Interestingly, simulation results bv lBurns et alJ (2008) indi¬ 
cate that non-cool core clusters (cool cores here are defined 
by a central temperature drop in the temperature profile) 
at low redshifts (z < 0.3) have less halos in their neigh¬ 
bourhood as compared to CC clusters, while the situation 
was exactly reversed at z > 1. Though burns et all (12008 1 
only discuss about clusters with virial masses greater than 
10 14 M 0 , the under dense regions of fossil systems relative 
to other systems would seem to qualitatively agree with the 
idea that the progenitors of fossils originated in a region 
with a higher density of halos that strongly contributed to 
mergers, and the final mass assembly currently resides in a 
relatively under-dense region. This is also quali tativ ely in 
agreem ent with the observational results of lHarrison et al.l 
( 2012 1 who base their conclusions via a study of the fossil 
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ellipticals and the simulation results of lDiaz-Gimenez et al.l 
( 2011 1. 

The above speculation however, still does not explain 
the ‘classical’ temperature profile seen for clll59p5531 and 
the high temperature fossils. One possibility is that the 
thermodynamic properties of the fossil systems are strongly 
dependent on the environment in which it was formed and 
it could be that these anomalous objects grew in a region 
relatively under-dense in halos that led to the formation 
o f a proper cool-co re via hierarchical, small-scale mergers 
((Burns et al.l I200S ) and led to relatively gentle AGN feed¬ 
back throughout. Secondly, we don’t rule out the possibil¬ 
ity of the high tem perature objects being falsely identified 
as fossils as the lEierenthaler fc Zeilingerl yfXhSf) sample is 
largely based on the SDSS catalogue and it has been known 
for some systems to be wrongly classified as fossils in simi¬ 
lar studies, on ly to be corrected late r w hen better data wa s 
available (e.g. IZarattini et al.ll20T3 vs. ISantos et al.ll2007h . 
C11159p5531 has a galaxy just outside the search radius 
w hich if included would no t satisfy the magnitude criterion 
of lYoevodkin et ah (1201 01 and be classified as a fossil in 
that study. Naively though, it seems unlikely that subtle 
changes in the magnitude gap would have a drastic impact 
on the ICM properties of these objects. 

4. Summary 

We have studied the core properties of 17 fossil systems 
from literature with data from the Chandra archives. The 
conclusions of our work can be summarised as follows: 

— Most fossil systems (82%) are identified as cool-core ob¬ 
jects based on at least two diagnostics. Interestingly 
though, there are indications that most fossil systems 
are WCC objects and that the population of non-SCC 
fossils outstrip SCC fossils. 

— Fossil systems show a large range of temperature profiles 
with many cool-core objects lacking the expected clear 
central temperature drop. 

— We analysed the X-ray emission coincident with the 
BCG for fossils with z < 0.1 and conclude that these are 
not the typical X-ray coronae reported for some other 
objects. We speculate that the stellar X-ray component 
(if present) has mixed with the ICM. 

— We performed a deprojection analysis using the SBPs 
for fossils with z < 0.05 and derived their thermody¬ 
namic properties. These objects lack group-sized cool- 
cores and show evidence for non-gravitational processes 
(AGN feedback in particular) based on the shallow pow- 
erlaw indices of their entropy profiles. 

— There are some indications that the normalisation of the 
Lx — T relation is higher for fossils than for non-fossils. 

— We speculate that early, dissipational major mergers led 
to a strong fuelling of the central SMBH leading to pow¬ 
erful AGN feedback and this could explain the lack of 
group-sized cool-cores in low-temperature fossil groups. 
There are however anomalous cases which lead us to 
conclude that this is probably an incomplete picture of 
the formation and evolution of fossil systems. 

In short, we have offered a glimpse into the nature of 
fossil systems through their core properties. In order to con¬ 
cretise these ideas, we will have to construct large, objec¬ 
tively selected fossils sample with high quality X-ray data; 


ideally with some mass/luminosity cut to distinguish be¬ 
tween fossil groups and clusters. A key point here seems 
to be to also improve our understanding of the weak cool- 
core cluster/group regime, which till date has been probed 
poorly. These observational results will also need backing 
from detailed, high-resolution simulations which factors in 
feedback processes that would help us unravel substantial 
features in the growth of fossil systems and the impact on 
the properties of the ICM. 
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Fig. A. 3. Temperature profiles of UGC 842 and RX J1331+1108. For RX J1331+1108, the blank-sky files might be 
under-estimating the background beyond rsoo- This should however not be a problem in the central regions. 
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Fig. A. 4. Temperature profiles of cll!59p5531 and cl!416p2315. 
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Fig. A. 5. Temperature profiles of RX J0825.9+0415 and RX J0801+3603. 












